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^\ ' We construct a field-theoretic description of two coupled spin-1 Heisenberg chains, starting with 

0^ . the known representation of a single spin-1 chain in terms of Majorana fermions (or Ising models). 

0\ • After reexamining the bosonization rules for two Ising models, taking particular care of order and 

disorder operators, we obtain a bosonic description of the spin-1 ladder. From renormalization-group 
and mean-field arguments, we conclude that, for a small interchain coupling, the spin-1 ladder is 
^ ' approximately described by three decoupled, two-frequency sine-Gordon models. We then predict 

I that, starting with decoupled chains, the spin gap decreases linearly with interchain coupling, both 

in the ferromagnetic and antiferromagnetic directions. Finally, we discuss the possibility of an 
incommensurate phase in the spin-1 zigzag chain. 
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I. INTRODUCTION 



Among the properties of spin ladders, the best known is the reduction of order as we go from a single spin-i 
Heisenberg chain to two coupled chains : the single spin-i chain is critical (its correlation length is infinite) whereas 
[ the spin-i ladder has finite-range correlations and an excitation gap, growing linearly with interchain coupling J±, 
^ ' at least for small J± (for a review and further references, see Ref. |l|). This may seem paradoxical because one would 
O ' naively expect that coupling two quasi-ordered chains would only increase the tendency to order, but a critical system 
, like the spin-i Heisenberg chain is easily sent off-criticality by a perturbation such as ladder coupling J±. In this 

paper we will study the corresponding spin-1 ladder (two coupled spin-1 chains), which is already disordered and has a 
T-H : finite gap at J± = 0. On the contrary, we will argue that the spin gap decreases as Jj_ increases from zero, and does so 
^ for both antiferromagnetic and ferromagnetic interchain couplings, thus giving the gap A{J±) a nonanalytic behavior 
(a cusp) at zero (cf. Fig. ^ below). We will arrive at this conclusion after obtaining a field-theoretic description of the 

■ spin-1 ladder in terms of six quantum Ising models, or alternately in terms of three boson fields. The motivation for 
QQ ] using bosonization is that it offers a safer description of the system at weak J±, valid for both positive and negative 

. J_L, and allows at the same time for a description of the spin-1 zigzag chain, in which frustration plays a role. Thus, 
^\ ' at small J±, this method is more general and reliable than a sigma model description. For a small antiferromagnetic 
0^ ■ interchain coupling, the drop in the gap as a function of J± was already ripticed in Monte Carlo simulations and 

' accounted for with a nonlinear sigma model description of the spin-1 ladder. □ 

I We will consider the spin-1 ladder as a perturbed critical model, so that the low-energy description of the system will 

■ be a perturbed conformal field theory. The critical model used as a starting point is a pair of decoupled biquadratic 
spin chains, with Hamiltonian 



Hq — |Sc(,i • Sc^i+i — (Sa,i ■ Sc(,i+i) I (1) 



] where Sa^i is a spin-1 operator at site i on chain cj-,J^ = 1,2). At this critical point the two chains are decoupled, 
■ each chain being described by an integrable modclElo which is equivalent in the continuum limit to a level-2 su(2) 
5—1 ' Wess-Zumino-Witten (WZW) modeLEl We then need to consider the following perturbation: 

a 

_ _ I 

Hi = {l + r^)J2 (So,. • S„^,+i)' + i Si.. • [(1 + + (1 - <5)S2,.+i] (2) 



I. a 



When (1 + 77) > 0, the first term brings us back to the Heisenberg point {rj = 0). The interchain interaction, of 
strength J±, is that of a ladder {S — ±1) or of a zigzag chain {S — 0). We will proceed by (i) constructing a continuum 
description of the interaction in term of WZW models and (ii) finding out the behavior of this perturbed WZW model 
by field-theoretic methods, mainly through representations in terms of Ising models (fermionization) and sine-Gordon 
models (bosonization). 
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FIG. 1. Schematic illustration of the coupled spin chains with the various couplings, normalized to the intra-chain coupling. 



This paper is organized as follows. In Sect. ||, we review the field-theoretic description of a sing le spin-1 chain, 
in particular its representation in terms of three Majorana fermions (or Ising models). In Sect. |ll, we write down 
a representation of two coupled spin-1 chains in terms of three bosons, using the bosonization formulas for pairs of 
Ising models given in Appendix ^ In Sect. the behavior of the spin gap as a function of interchain coupling 
is inferred from this bosonized description. In Sect. 0, the spin-1 zigzag chain is considered instead, and a weak 
interchain coupling is argued to cause a short-distance incommensurability, i.e., a displacement of the minimum of 
the one-magnon spectrum from q — it. 



II. CONTINUUM DESCRIPTION OF THE SINGLE SPIN-1 CHAIN 



A. Phase diagram 

Let us first review the phase diagram of the biquadratic spin-1 chain:!! 

H = Y,{^^■^^+l+V{^^■^^+lf} (3) 

i 

At 77 = —1, the Hamiltonian is integrable and has gapless modes at fc = and k — tt. It is also integrable at 77 = 1 
and has then gapless modes at fc = and k = ±27r/3. If 77 < —1 we have a dimerized phase characterized by two 
degenerate ground states with a finite gap. On the other hand, in the interval 77 € (—1, 1) the spectrum has a singlet 
ground state with a finite, gap. This is the so-called Haldane phase, characterized by the spontaneous breakdown 
of a Z2 X Z2 symmetry.LrB This breakdown implies a four-fold degenerate ground state in an open chain, but these 
different ground states differ only by the spins at the ends of the chain, and in this sense they are equivalent in the 
thermodynamic limit. The excitations are solitons switching from one ground state at a; ^ — c» to another ground 
state at a; — > 00. Related to the symmetry breaking is a dilute antiferromagnetic order; schematically: 

+ 0---0-0---0 + 0---0- (4) 

This order is defined by an alternation of sites with 3.^ = 1 and Sz — —1, with some 5'^ = sites in between. It can 
be measured by the so-called string order parameter:!!! 

O^^ lim (s:,expl^7^J2sAs^) (5) 

\ \ k=n+l / I 

This order parameter and the gap are maximal at 77 = 1/3, where the VBS-like ground state is exactly knownB The 
gap grows monotonically from 77 = — ltor7 = without phase transition, and thus we may consider the Heisenberg 
point (?7 = 0) as a perturbation of the critical point (77 — —1). 

Note also that incommensurability develops starting at 77 w 0.4: the peak in the spin-spin correlation function 
moving from fc = tt to /c = 27r/3 at n-i= l.!l3 This last transition point is described by an SU(3) generalization of the 
Kosterlitz-Thouless phase transition.lii! 



B. Field-theoretic description 



The critical point (77 — —1) is equivalent, in the low-energy limit, to a conformal field theory: the su(2) Wess- 
Zumino-Witten (WZW) model at level /c = 2, plus a marginally irrelevant perturbation.B This WZW model contains 
two scaling fields : a spin doublet gmn {m^n € {—5, 5}) with left and right conformal dimensions (^, ^) and a spin 
triplet (m, n G {—1,0,1}) with dimensions (5,5). They are respectively 2x2 and 3x3 matrix fields. The 
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link between thej^pin chain and the WZW model is given by the following representation of the spin operators in the 
continuum limit 

-S^{x) = i {rix) + r{x)) + i-ir/^oQg'^ix) (6) 

where ao is the lattice constant, Q a nonuniversal constant, J° and J° are the right and left su(2) currents and g°' is 
defined in terms of Pauli matrices as 

5" = -^Tr {a^g) = E <nn9nm (7) 

The currents (J", J") and the field correspond to the soft modes of the spin chain near fc = and tt, respectively. 
For 1 + 7] not too large, the spin chain may be described by the above WZW model, plus the following perturbation:tl 

£1 = mTr$ - AiJV" (8) 

where a summation over repeated indices is implicit and —m is proportional to (1 + 77) (m is negative in the Haldane 
phase). The second term is the marginally irrelevant perturbation alluded to above (if Ai > 0). On the other hand, 
the first term (Tr$) is relevant, with scaling dimension 1, and leads to a gap proportional to \m\ oc |1 + ri\. j— . 

There is an interesting equivalence between the k = 2 su(2) WZW model and three quantum Ising modcls,E3 and 
so we will not have to deal with the WZW model directly. This equivalence is defined by the following relations: 

— 2 — — i — — 

J"^ = —^eabclpb'ipc = —^eabcipbi>c (9) 



$1 = — {-^pi+i-ip2) 



^0 = CV'S 



7!'* 



#2) 



(10) 



$0 = CV's 



$„i = (V'l - #2) 



(11) 



= V2CTiCr2Cr3 g" = -V2aa^J.a+l^J■a+2 (12) 

where the latin index goes from 1 to 3; ipa and ipa are respectively the right and left fermions associated with each 
Ising model (see Appendix ^). o-q and Ha are the order and disorder fields of each Ising model. The 3x3 matrix field 
^nm is here factorized as = The constants ( and ( are such that their product is CC = Note that our 

relations differ slightly from those given by Fateev and Zamolodchikov.E^I The action of the WZW model in imaginary 
time becomes simply that of free Majorana fermions: 

S'wzw = ^ J dxdr {ipadi'a + "^adi^a) (13) 

where d — {dr — idx)/2 and d = [d^ + idx)/2 (in order to lighten the notation, the characteristic velocity v of the 
WZW model has been set to unity). The perturbation (^) becomes: 

Cl = m-)pa4'a - Xl1pa1pa'>Ph1pb (14) 

Except for the marginally irrelevant term, the spin chain is thus equivalent to three Majorana fermions of mass rpu 
This description of the spin-1 chain has been used to study the effect of a magnetic field on the low energy spectrum.E^ 
The representation (||^^2|) of the WZW fields is invariant under the following changes (for a = 1,2,3 simultaneously): 

Ipa -Ipa Ipa -'4'a A^a ^ "Ma ^ Ca (15) 

This is related to the absence of fermionic field in the WZW model. This 'gauge' symmetry accounts for the expected 
degeneracy of the ground state near the critical point in open chains. Specifically, recall that m < in the Haldane 
phase. In our formulation, this corresponds to the disordered phase of the Ising models (see Appendix ^ and the 
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expectation value of the disordered operators is nonzero: (/ia) 7^ 0. In this phase each Ising model has a doubly- 
degenerate ground state, associated to different spin configurations at the ends of the open Ising chain. The two 
ground states differ in the sign of (/in). For the spin chain, this degeneracy implies an apparent eight-fold (8 — 2^) 
degeneracy, but the gauge invariance ( [l5| ) reduces this to a physical four-fold degeneracy. These different ground states 
come from the breakdown of the hidden Z2 x Z2 nonlocal symmetry alluded to above, and are physically equivalent 
in the thermodynamic limit. In this Ising model description of the spin chain, the elementary excitations are-|kinks 
switching from one value of (/ia) at a; ^ —00 to its opposite at x — > cxd. On the other hand, Fath and Solyomcl have 
shown that the excitations of the Heisenberg model are solitons connecting the ground states with different values of 
the string parameter . We are thus led to identify these solitons with the kinks of the Ising model. 

We can do the same exercise for to > (or ry < —1). We are now in the ordered phase of the Ising models: 
(fa) ^ 0. Such an expectation value is already invariant under the gauge change ([T^) and therefore there are really 
8 physically different ground states for the open chain. A hidden Z2 x Z2 symmetry breaking is again expected and 
so these 8 different ground states will be locally equivalent to two distinct ground states in the thermodynamic limit, 
corresponding to the expected dimerized state. 



III. BOSONIZATION 

Using the continuum description (^) of the spin operators, we obtain the following Lagrangian density from the 
Hamiltonian (|l|,§), in terms of WZW fields: 

C = /:wzw[<7] + -Cwzwb'] + mTr$ + TOTr$' - Ai(JV° + J'" J'") 

-f A2(JV° + rJ'"-) + XsirJ"" + JV") + Ag^^"^ + p (5°9^<?'° - (a:,5°)<?'°) (16) 

The unprimed fields correspond to the first chain and the primed fields to second chain. The first three terms represent 
the intrachain interaction and the last four terms the interchain coupling. The interchain couplings A2, A3, A and p 
are respectively proportional to J_l, J_l, J±S and J±{1 — S^) at high-energy, but they rennrmalize differently towards 
low-energy. The last term has been omitted in previous work on the zigzag spin chain,t£Hl3 but since it respects the 
symmetry of the microscopic model, it is not forbidden and we expect p to be nonzero. Note that p must be zero for 
the spin ladder, whereas A vanishes for the pure zigzag chain. In the following we will consider the Haldane phase 
only so that to is negative. 

The Lagrangian ( |l6| ) is difficult to study in terms of WZW fields. The simplest information we may extract from 
it is the scaling dimension of the various perturbations, from those of the various WZW fields. Thus, the interchain 
couplings A2, A3, A and p respectively have scaling dimension 2, 2, |, and |. Moreover, the couplings A2 and p have 
conformal spin. By itself, a relevant coupling g of scaling dimension 7 < 2 and zero conformal spin is expected to 
produce a gap of order A ^ g^^^^^'^\ Thus, at the in-chain critical point (77 = —1), the interchain coupling A would 
open a gap of order 

A(A) - A^/^ (17) 

in the spin-1 ladder. 

However, far from the critical point, the WZW model is of little help in predicting the behavior of the gap and the 
fermionic language seems more appropriate. Using the representation (p|-p^), we can express the Lagrangian density 
( |l6| ) in terms of Majorana fermions, order and disorder fields. Unfortunately, the resulting expression is not easy to 
study since it contains a mixture of fields that are mutually nonlocal (the order and disorder operators). 

An interesting way to deal with the Lagrangian ( [l^ ) is bosonization. The 2D Ising models may be bosonized by 
pairing them (see Appendix ^). The natural way to bosonizc the ladder is to pair an Ising model describing one 
chain with its twin on the other chain. Using the relations (p| , ( [l6| ) and (All), we obtain the following Lagrangian 
density for two coupled spin-1 chains: 

C = Co + Ci+C2 + C3 + Cp + Cx 

£0 = X! {l^:^ [{dr<Paf + {d.j,ipaf] " 2to COS 
a=l,2,3 

£1 = I6A1 ^ ( COS ipa+1 COS Lpa+2 + COS Oa+l COS Oa+2) 
a=l,2,3 

£2 = 4A2 ^ {dr^Pa+ldr^a+2 ~ dxiPa+ldx^a+2) 
a=l,2,3 
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£3 = -8A3 ( 

sin ipa+i sin ipa+2 + sin 9a+i sin 9a+2) 

a=l,2,3 

C\ = 4v2A 2^ cos — sm sm 

a=l,2,3 



£p = -4V2p ^ cos6ia 



fa . fa+l . fa+2 . Va . ( Va+1 + 

cos — — sm — - — sm — sm — — sm 



2 2 2 2 V 2 



(18) 



a=l,2,3 

where is the boson dual to ^pa- To shorten the expression, we have adopted a periodic condition on th e index a, i.e., 



a + 3 = a. The twist term £p, the trickiest to bosonized, has been inferred from the representation (A12) of the stress- 
energy tensor for each Ising model, plus the usual operator product expansion (OPE) between the energy-momentum 
tensor and a conformal field. 

Thus, we have transformed the problem into a system of three perturbed sine-Gordon models, although the simul- 
taneous presence of the bosons kp^ and of their dual fields da makes some perturbations nonlocal. However, as we 
will see, the most relevant perturbation is local and makes the problem tractable in this language. Note that our 
normalization is such that cos(/3(^a) is marginal for /3 = v^, and thus bound states appear in the sine-Gordon model 
for /3 < 1. Also, for (i < •\/2, the ipa Pa + '^n symmetry is spontaneously broken and we have to consider fiuctuations 
around one of the minima of the potential. However, it is important to keep in mind that our bosonization procedure 
is from the start invariant under the translation ipa lyS^ + 47r, and this 47r-periodicity property must be regarded as a 
constitutive constraint imposed on the sine-Gordon models. Thus, each sine-Gordon model in Co has two inequivalent 
ground states, associated with the minima ipa = ±7r of the potential (for m < 0). The spontaneous breakdown of the 
symmetry ipa — > 1^0 + 27r implies a nonzero expectation value for the operators siniy9a/2 (the disorder operators fia) and 
cosiy9a. Moreover, this breakdown becomes explicit if the perturbation Cx or Cp is added. This symmetry break ing o f 
the three sine-Gordon models corresponds in fact to the hidden symmetry breaking in the spin-1 chain (cf. Sect, A| ). 
Therefore the different choices of the ground state {cpa — ±7r) are equivalent, since the different ground states of the 
spin-1 chain are equivalent in the thermodynamic limit. Finally, let us recall that the elementary excitations of each 
sine-Gordon model have finite mass and correspond to the kink and antikink connecting the two different ground 
states. The charge conjugation changing kink into antikink corresponds to the following transformation: 

1^9(2;) 2tt — (p{x) mod 4tt (19) 

The presence of nonzero expectation values for the operators sin(/?a/2 and cos ipa implies that more relevant terms 
may be generated from the perturbations ([l^. If we call a the expectation value of sm.ipa/2, say in the (p — tt ground 
state, and ai that of cosipa^ then 

COS — sm sm = ±a cos h ffuctuations 20) 

2 2 2 2 ^ ^ 

The sign depend on the relative choice of the grpuiid state for Pa+i and Pa+2- The expectation values a and ai have 
been obtained by Lukyanov and ZamolodchikovEZlEa and are proportional to ItoI^/"* and m respectively. Keeping only 
the most relevant terms and neglecting the fluctuations of sin 1^9/2 and cos 1^9 around these expectation values, we find 
the following effective Lagrangian: 

^oS=J2{^ [i^rPaf + [d^cPa?] ~ (2m - 16Aiai) cos Pa ± ^^\o? COS 
a 

T'iV2pacos9a cos ^ - cos - cos } (21) 

At this level of approximation, we have three perturbed sine-Gordon models - mutually coupled only if p 7^ - 
and the sign of the interchain coupling can be incorporated in the choice of ground state. Thus, a ferromagnetic or 
antiferromagnetic interchain coupling would have the same effect. Note that the couplings A and p break the charge 
conjugation symmetry ([l9|). 



IV. BEHAVIOR OF THE GAP IN THE SPIN LADDER 



Let us first consider the spin-1 ladder, which corresponds to p 
three decoupled, two-frequency sine-Gordon models: 



0. The effective Lagrangian (21) then reduces to 
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FIG. 2. Evolution of the sine-Gordon potential cos if + Acos((^/2) for A = 0, 0.1 and 0.5. The mass of the lowest-energy 
kink decreases linearly with A. 



Aad [(dr^a)^ + {d:.Vaf] - M COS (^a ± A COS ^ | (22) 

where M = (2m — 16Aiai) and A = iy^Xa^. This Lagrangian has been studied by Delfino and Mussardo,lli who 
conclude that, as A increases from zero, one of the kinks from each model becomes more massive, whereas the other 
one becomes less massive. This is easily understood by considering the evolution of the potential V{(pa) as a function 
of A (cf Fig. 1^): the soliton having to bridge the potential barrier from ip ^ tt to ip ^ 3tt (towards the right) has 
a lower energy than the soliton going from tt to (/s ~ — tt = Stt (towards the left). Which kink sees its mass 
decrease depends on the sig n of t he perturbation, but the net result is the same whatever this sign is. With the help 
of sine-Gordon forin_iactors,ll2H13 we can ascertain how the kink mass varies with A. At first order, the variation of 
the mass squared isxB 

Sml « \A\Faa{i7r) (23) 

where the form factor F is 

FaM = {0\sm^\a{f^,)a{7j2)) (24) 

where a and a represent the kink and antikink and 771^2 are the associated rapidities {rj — rji — 772)- From Ref. we 
extract the following expression: 

Fa-airj) = - (e*^/2) e"/2 - (e-'W2) ^-^/^ 

la) 21/6^3 e-i/42 cosh f (25) 



where is the mass of the kink and A « 1.282427 is the Glaisher constant. From this result, we see that Sml 
vanishes at first order. We thus expect it to be proportional to A^, This is compatible with the semiclassical result 
that the variation of the mass of the kink is proportional to the variation of the height of the potential. We thus 
conclude that 

Sma (X A (26) 



According to the analysis of Ref. 19, a single two-frequency sine-Gordon model has a c = ^ Ising fixed point, in 
addition to the Gaussian fixed point at Af = and A = 0. Since the scaling dimensions of M and A at the Gaussian 
fixed point are respectively 2 — 1 = 1 and 2 — i = |, the ratio C = K/M''/^ is invariant under RG flow and is in fact 
a control parameter which tells us how far we are from the Ising fixed point, characterized by a critical value Cc. At 
this value, i.e., at A = QcM''^'^, the light kinks have exactly zero mass. If we return to an Ising-model description of 
the system, we can understand intuitively how this flow happens: The effective Lagrangian (E3) corresponds to six 
2D Ising models coupled pairwise by the following interaction: 



gap A 




J± 



FIG. 3. Conjectured dependence of the spin gap A upon the interchain coupling J± in the spin-1 ladder. 

Thus, the excitation such that a{x) is parallel to (j'{x) will have a lower mass if A > (a similar reasoning holds when 
A < 0, by changing the sign of <t'). When A is large enough, <t must be parallel to a' and this parallel configuration 
defines a new Ising model, whose critical point occurs at some value of the ratio A/M^/'*. 

Whereas the mass gap vanishes at a finite value of A in the effective Lagrangian (^2|), this does not occur in 
the spin-1 ladder, because of the marginal terms £2 and £3 and of the (neglected) fluctuations of Eq. (po|). We 
nevertheless expect that the spin-1 ladder should flow close to this fixed point and that the gap should decrease 
linearly with a weak interchain coupling, both on the ferromagnetic and antiferromagnetic sides (with the same 
slope). This is to be compared with the Monte-Carlo data of Fig. 3 of Ref. H, which illustrates this drop in the gap, 
for an antiferromagnetic interchain coupling only. That the gap drops on the ferromagnetic side is not surprising, 
considering that (i) the ladder becomes equivalent to a spin-2 chain at large ferromagnetic coupling and (ii) the gap of 
an antiferromagnetic Heisenberg chain with integer spin s decreases with s. On the antiferromagnetic side, the drop 
in the gap may be roughly understood as a greater tendency towards magnetic ordering, and should be maximum 
when J_L ~ 1. Beyond J_l ~ 1, the gap must eventually turn around and increase linearly at large J±, since the 
lowest-lying excitations are then rung triplets, costing an energy J^. The gap A{J±) is then conjectured to have a 
cusp-like maximum at J± = 0, a, peculiar nonanalytic feature, as illustrated schematically on Fig. 



V. THE ZIGZAG SPIN CHAIN 



The zigzag spin-1 chain corresponds to d — 0, and thus A = 0, /? 7^ 0. The effective Lagrangian is then 

-Czigzag = X! 1 ^ [(<9t¥'q)^ + {dxf)'^] - (2m - 16AiQ;i)cOS((9a 

=F4V2pa cos {0a) 



fa ifa+l fa+2 

COS COS cos 

2 2 2 



(28) 



This Lagrangian is not easily analyzed. Let us go back to the fermionic representation of the twist term by order and 
disorder fields: 



£p = 2p^ 

(^a^J'a+l^J■a+2^xicra^J■'a+l^^'a+2) 

a 

With {fiaiJ,'a) — —V^ia (cf Eq. All), the most relevant term will be 

Cp w -Aa'^paad^a'a 



(29) 



(30) 



We will now study the effect of this approximate representation of the twist term by considering the corresponding 
lattice model (see Appendix ^). Let us map the order fields in the following way: 

aa{x)^a:{n) a',{x) (n + (31) 
With the representation ( ^ ) for the twist term, the system is described by the following Hamiltonian: 

^ = E {-'^a(^V2) - «<(n/2)<(n/2 + 1)} - W ^ K{n + 1/2) - <(n - 1/2)] (32) 
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where k is related to the constant m (i.e. —1 — rj) by the relation k = 1 + aom, where qq is the lattice constant. Thus 
K = 1 for m = (ji = —1) and tends to when i] grows. To bring this Hamiltonian to a more familiar form, we 
perform a rotation of tt around the z axis of the spin operator at every other site, on each chain. This changes the 
sign of (T^ on those sites and gives the Hamiltonian a slightly different form: 

H = Y. + ^<{n/2)a:{n/2 + 1)} + 4^^^ ^<(n) K(n + 1/2) + a^,{n 1/2)] (33) 

n,a n,a 

The Hamiltonian defines the quantum ANNNI model. Together with its two-dimensional, classiGaJ, 

counterpart (cf Refs EC^g^L it has been extensively studied by a pSiariety of methods: mean, jfield theory,c^ 
Monte Carlo srapkytiDns^OLlI Muller-Hartmann-Zit|tartz approximation^ perturbative expansions ,E^Ej free fermion 
approximationEjlEJ^Ell and exact diagonalizations.EircS The phase diagram for the classical model is shown in Fig. |^. 
In the scaling limit, the temperature T of the classical model model is related to the mass to by T = (1 — aom)Tc = 
(2 — k)Tc. The NN coupling is proportionnal to the interchain coupling p and J2 to n. Thus, the case of small zigzag 
interaction corresponds to the limit of small Ji. The different phases are the following: ferromagnetic (F), para- 
magnetic commensurate (PC), paramagnetic incommensurate (PI), incommensurate critical phase (IC, also called 

"floating phase") and antiphase (A) of alternating pairs (+ H h + • • •). A disorder line found by Peschel and 

Emer5iH divides the PC and the PI phase. 

We conclude from this phase diagram that incommensurability will arise in the spin-1 zigzag chain as soon as the 
interchain coupling is nonzero (the model (^3|) is then in the far-right of the PI phase) . One premise for this deduction 
is that the incommensurabilityof the Ising spins ((Tq) is reflected in the correlation of the spins of the quantum chain; 
this comes from the relation ( p^ . Note that increasing p brings us from infinity on the phase diagram ^ towards the 
origin, along a straight line. One could expect such a line to go through other phases (like the IC phase) at some 
point. However, we should note that the omission of the fluctuation of Ha in Eq. (|3^) is valid only when T — Tc is 
large compare to p. 

Moreover, we can have an idea of how the inconuaensurability develops as a function of p. A recent analysis 
using a high-temperature expansion and bosonisationEJ shows that in the limit of very strong next-nearest-neighbor 
interaction in the ANNNI model, the incommensurability is proportional to p/k. Explicitely, in the high temperature 
limit, the incommensurate wavevector is given by: 

,. = ±HfMi±^±4±:_i (34, 

2k(1 — K) 

where the dots stand for higher power of k. The ± sign are respectively associated to the correlation function of the 
combination ±a^{n) + a^{n^ 1/2)- This result that—the incommensurability is linear with the interchain coupling 
confirms the one obtained by a semiclassical analysis.E3 



ACKNOWLED GMENTS 



D.A. thanks Philippe Lecheminant for useful discussion. This work was partially supported by NSERC (Canada) 
and by le Fonds FCAR (Quebec). 



8 



APPENDIX A: THE 2D ISING MODEL 



In this appendix, we review briefly the correspondence of the Ising model with fcrmions, the conformal structure 
of the model and we indicate a set of careful bosonization formulas for a pair of Ising models. 



1. Definitions 



As is well known, the 2D statistical Ising model is equivalent to a quantum Ising chain in a transverse field, with 
Hamiltonian 



where tr^'^''^ are the Pauli matrices. The Hamiltonian (Al) can be diagonalized through a Jordan- Wigner transfor- 



mation followed by a Bogolubov-Valatin transformation. The solution shows that (erf) 7^ if A < 1 and (erf) = 
otherwise. Thus A = 1 is the critical point. A peculiarity of this model is the existence of a duality transformation 
mapping the ordered phase to the disordered phase and vice-versa. Under this transformation the spin operators cr" 
are mapped to the so-called disorder operators, defined on links (dual lattice) by the following relations: 

M-+i/2-n^i (A2) 

Let us apply the Jordan- Wigner transformation on the dual lattice. The fermion creation and annihilation operators 
are defined as: 

Cj+1/2 = ^7+1/2 exp Y ^('^fc+i/s - 1) 



4+1/2 = A^j>i/2 exp E(^fe+i/2 - 1) (^3) 



-ZTT 

,± 



where [i = (^^ ± z^^)/2. The fcrmions correspond to the kinks in the original formulation. Indeed the fermion 
number on link j -|- 1/2 is 

•=^+1/2^^+1/2 - 2 ' ^ ^ 

i.e., there is no fermion on the link if the spins on j and j + I are parallel and one if they are antiparallel. Note that 
the order parameter has a bosonic character, whereas the disorder parameter fj,^ is fermionic. This is easily seen 
from the following equivalence: 

Aij + 1/2 = 4+1/2 -*7r^cj'_^i/2Cfe+i/2 +Cj + i/2exp 

\ k<'J j \ k<] 

<7j ^(tIn exp ( iiTT J2 4+l/2'^fe+l/2 I (A5) 



2. Continuum limit 



The critical point of the 2D Ising model is equivalent, in the continuum limit, to a free, massless Majorana fermion: 
a conformal field theory with central charge c = 1/2 and three conformal families: the identity operator, the energy 
operator e, and the the spin density (or order) operator a. The use of complex coordinates z — t + ix and z = t — ix 
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is standard, along with the complex derivatives d = dz = (9r — idx)/'2 and d = dg = {dr + idx)/^ (the notation 
used is that of Ref. ^6|) . The energy density operator may be expressed in terms of the chiral components of the 
Majorana fermion as e = ii/jip. The order field a is the continuum limit of the spin operator erf, and a fermionic 
disorder field /i may be introduced as the continuum limit of the disorder operator Mi+i/2- The field /i has the same 
scaling properties as the field a, but is nonlocal with respect to a. The conformal transformations are generated by 
the energy-momentum tensor, whose chiral components are T — and T = — i^/jQ-f/). All these fields have the 

following short-distance products, or operator product expansion (OPE): 



<t{z, z)(j{w, w 
fi{z, z)^{w, w 

<7{z, z)^{w, W 



1 



z — w 

1 

z — w 

1 



+ 2(z- w)T(w) 
+ 2{z-w)T{w) 



— + -\z-w\"\{w,w) 

1 ^171 - \\z-w\^'^e{w,w) 

\z~w\^i^ L 

7(z — w^l'^^iw) -f 7*(£ — wYl'^^{w) 



V2|z-w|i/4 



/i(z, z)a{w 


w) ' 


7*(z - w)^/^^(u)) -1- 7(z - 


V2|z-u;|i/4 


ip{z)a{w 


w) ' 


7 

V2iz~wy/-^ ' 


ip{z)fi{w 


w) ' 


7* 

^ —;= t(^{w,'w) 

V2{z-wY/^ ^ ' 


V'(z)cr(w, 


w) ' 


^V2(z--^)V^^^"'"^ 


■0(z)/i(w, 


w) ^ 


—aiw^w) 

^{z-wy/^ ^ ' 



^1/2, 



where 7 = expi7r/4 (or, equivalently, exp— z7r/4). 



3. Bosonization 



(A6) 



Two Ising models form ac = 1/2-1-1/2 = 1 conformal theory. We therefore expect a representation of the different 
fields in terms of a free boson defined by the action 



5 = — / d-xdr 
Svr 



(A7) 



Our choice of normalization (I/Stt) simplifies the exponentials and circular functions appearing in the sine-Gordon 
theory. The OPE of the boson field and of its (normal-ordered) exponentials are 



d(p{z)dip{w) 



iaLp(z) i/3ip(iu) 



(z — wY 



{z-wT'' e 



ap i{a+f3)(p{w) 



The boson field can be separated-Luto chiral components: (p{x, t) 
expansion in radial quantization£3 



(A8) 

(z) . These fields have the following mode 



(z) = q-ip\ogz + i^^ 



{z) = q — iplogz + 



k^Q 



dkZ 



auz 



(A9) 
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where the operators p, q, an satisfy the following commutation relations: 

[<l,p] = * [an, am] = n5n,m 

with similar relations for the left-moving (barred) operators. j— , 
A faithful representation of the Ising fields is then given by the following relationsiEZi 



(AlO) 



ajj,' 
/iff' 



1 

71' 



'V2' 
i 

71' 
1 

71' 



1 

71' 

-10 



71" 



-ITTp 



71 



-nrp 



71 



V2 COS ^ 



-iv^ sin — 

■ /7T & 
'tV2 cos - 



= -V2: 



I sm ■ 



(All) 



where 6 ^ (j) ^ cj) + 2'np is the field dual to Lp (the operator 2'Kp is added to ensure proper anticommuation properties). 
This representation leads to the correct OPE ( |A6| ) between the Ising fields. The phase factor e*'^^ is similar to the 
phase factor in the Jordan- Wigner transformation, p being to the number of left fermions. Only its odd or even 
character matters. We note the natural periodicity property Lp ^ Lp + Att and 9 ^ 9 + An oi this representation. 

The energy-momentum tensors T and T' of the two Ising models, along with their antiholomorphic counterparts, 
are bosonized as follows: 



T(2) - T'(z) = 4\/2e2"Pcos (2(/)(z)) 
f{z) - f'{z) = -4\/2e2^*Pcos(2<^(z)) 



(A12) 



This relation is useful when bosonizing the twist term (the last term of Eq. ( [16| ) 
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